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ABSTRACT 

We present measurements of projected autocorrelation functions Wp{rp) for the stellar 
mass of galaxies and for their light in the U, B and V bands, using data from the third 
data release of the DEEP2 Galaxy Redshift Survey and the final data release of the 
Sloan Digital Sky Survey (SDSS). We investigate the clustering bias of stellar mass 
and light by comparing these to projected autocorrelations of dark matter estimated 
from the Millennium Simulations (MS) at z = 1 and 0.07, the median redshifts of our 
galaxy samples. All of the autocorrelation and bias functions show systematic trends 
with spatial scale and waveband which are impressively similar at the two redshifts. 
This shows that the well-established environmental dependence of stellar populations 
in the local Universe is already i n pla ce at z = 1. The recent MS-based galaxy for- 
mation simulation of iGuo et al.1 (|201lD reproduces the scale-dependent clustering of 
luminosity to an accuracy better than 30% in all bands and at both redshifts, but 
substantially overpredicts mass autocorrelations at separations below about 2 Mpc. 
Further comparison of the shapes of our stellar mass bias functions with those pre- 
dicted by the model suggests that both the SDSS and DEEP2 data prefer a fluctuation 
amplitude of cts ^ 0.8 rather than the cts = 0.9 assumed by the MS. 

Key words: galaxies: clusters: general ~ galaxies: distances and redshifts - cosmology: 
theory - dark matter - large-scale structure of Universe. 



1 INTRODUCTION 

The two-point correlation function (2PCF) has long served 
as the primary way of quantifying the spatial distribution of 
galaxies in our Universe. Measurements of 2PCF for differ- 
ent classes of galaxies in the local Universe have been car- 
ried out with high accuracy thanks to the large redshift sur- 
veys assembled in recent years, in particular the two-degree 
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field galaxy redshift survey (2dFGRS; IColless et aD 120011 ) 



and the Sloan Digital Sky Survey fSDSS: lYork et al.|]2000l i 



These studies have established that 2PCFs depend on a va- 
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IWang et al]|2010l : IRoss et al.ll201l] : IZehavi et al.ll201ll ). Re- 
cently, there have also been studies of galaxy clustering 
at higher redshifts, which are usually based on deep sur- 
veys (up to z ~ 1 — 1.5) covering small areas (<1.5 deg^) , 
e.g. the DEEP2 Galaxy Red shift Survey (jPavis et al.ll2003l : 
ICoil et al.l200 4a'.'2006'. 2008V the VIMOS-VL T Deep Survey 
dLe Fevre eTal. 2005 : Polio ct al, 20Q5 , 20061; iMeneux et al.l 
20061. l2008l)~and the zCQSMOS Survey (iLillv et all l200l 



2003 : iMeneux etahl |2009| : Ide la Torre et all bOllaf T Mea 



surements of 2PCFs for different classes of galaxies and at 
different redshifts have provided powerful quantitative con- 
straints on models of galaxy forma t ion and evolution (see for 
exarnnlelPe Lucia fc BlaizotllioOTi : iLi et al.|[2007l : IGuo et all 
l201ll : ]de la Torre et al.ll2011bl ). 

A number of recent studies have further investi- 
gated the dependence of galaxy clustering on physical 
properties by measuring weighted or marked 2PCFs (e.g. 
Beisbart fc Kerschei I2OOOI : iFaltenbacher" et al.l l2002l: IShethI 



2005':ISkibba ct al.'2006':'Sk ibba fc Shethl2009l : lskibba et al 

200 9; Li & White 2009, 20i3)- Simply speaking, this statis- 
tic is estimated using exactly the same methodology as the 
one used for the traditional 2PCF, except that each galaxy 
in the real sample and/or in the random sample is weighted 
by one of its physical properties such as stellar mass or 
luminosity at a given band. When compared to the tra- 
ditional 2PCF, this alternative two-point statistic has the 
advantage that it makes use of the whole galaxy sample 
and thus minimizes the sampling and large-scale structure 
noises. As found in lLi fc White! ( 20091 ). a particular virtue 
of the weighted correlation functions is that the correlation 
signals are dominated by contributions from galaxies in a rel- 
atively narrow mass range, around the characteristic mass 
of the stellar mass function. Thus, the statistic is robust to 
incompleteness of the observed sample at the two mass ex- 
tremes. In addition, the weighting scheme used for estimat- 
ing the statistic makes it a direct measure of the clustering 
of stars on scales larger than those of individual galaxies, 
rather than the clustering of the host galaxies as probed by 
the traditional 2PCF. This statistic thus provides a compact 
and accurate way to characterize the distribution of stellar 
populations over large ranges in spatial scale. 

Using a sample of almost half a million galaxies from the 
Sloan Digital Sky S urvey (SDSS). iLi fc White, (,20091 . here- 
after Paper I) and iLi fc Whit j (|2010l . hereafter Paper II) 
estimated projected autocorrelation functions Wp{rp) both 
for the stellar mass of galaxies and for their light in the 
five SDSS photometric bands. All of the autocorrelation es- 
timates are extremely well described by power laws over 
the full non-linear range 10/i~^kpc< rp < 10/i~^Mpc. Lu- 
minosity is found to cluster less strongly than stellar mass 
in all bands and on all scales. The autocorrelation function 
of luminosity varies systematically with wavelength in both 
amplitude and slope, with the reddest band (the 2 band) 
showing the highest amplitude and the steepest slope, indi- 
cating that the z-band light is the closest proxy for stellar 
mass in terms of clustering properties. These trends pro- 
vide a precise characterization of the well-known dependence 
of stellar populations on environment. In combination with 
autocorrelation functions of luminosity- and stellar mass- 
selected galaxy samples, as well as accurate luminosity and 
stellar mass functions, these results provide tight constraints 
on galaxy formation, which are a major challenge for cur- 



rent models (see, for exa mple, the significant discrepancies 
in recent work bv lGuo et al. 2011, based on the Millennium 
Simulations). 

In this paper, we extend the work of Papers I and II to 
higher redshifts ( z ~ 1) using data from the DEEP2 Galaxy 
Redshift Survey (jDavis et al.ll2003l ). We use the methodol- 
ogy of Papers I and II to compute projected stellar mass 
autocorrelation functions for DEEP2 galaxies, as well as 
projected autocorrelations for their luminosity in the rest- 
frame U and B bands. In order to compare the two sur- 
veys directly, we re-analysize the SDSS, computing lumi- 
nosity autocorrelations in the U, B and V bands, where 
luminosities in these bands are estimated from SDSS data 
through a Bayesian technique based o n a spectral energy 
distrib ution library constructed from the lBruzual fc Chariot! 
(|2003! . hereafter BC03) population synthesis code. We com- 
pare these results with the autocorrelations of dark matter 
at z = 0.07 and z — 1, the median redshifts of the SDSS 
and DEEP2 galaxy samples used here, in order to under- 
stand the bi as of stellar mass and light. By comparing with 
the model of !Guo et al.! (|201ll ). we investigate quantitatively 
how well current treatments of galaxy formation reproduce 
the clustering evolution of stellar mass and light. Finally, we 
discuss the possibility that the measured shape of the mass 
autocorrelation functions can be used to estimate the value 
of the mass fluctuation amplitude parameter as- 



2 DATA 

2.1 SDSS galaxy sample 

The low-redshift galaxy sample used in this study is 
a magnitude-limited sample constructed from the final 
data release (D R7; lAbazaiian et al.! 12009! ) of the SDSS 
(|York et al.ll2000l ) and is exactly the same as that used in Pa- 
pers I and II. This sample contains 482,755 galaxies located 
in the main contiguous area of the survey in the northern 
Galactic cap, with r < 17.6, —24 < Mo.i^ < —16 and spec- 
troscopically measured redshifts in the range 0.001 < z < 
0.5. Here r is the r-band Petrosian apparent magnitude, 
corrected for Galactic extinction, and Afo.i^ is the r-band 
Petrosian absolute magnitude, corrected for evolution and 
_ft'-corrected to its value at z = 0.1. The apparent mag- 
nitude limit is chosen in order to select a sample that is 
uniform and complet e over the entire area of the survey (see 
iTegmark et al.![2004! ). The median redshift of this sample 
is z = 0.088, with 10% of the galaxies below z=0.033 and 
10% above z=0.16. As shown in Paper I (see their fig. 4) 
the autocorrelation of stellar mass is dominated by contri- 
butions from a narrower and slightly lower redshift range 
(with 10% of the galaxies below z — 0.025 and 10% above 
z — 0.12), with a median redshift of z = 0.067. We thus 
take 2 = 0.067 as the effective median redshift of this sam- 
ple when comparing with dark matter autocorrelations and 
semi-analytical model predictions. 

We use a Bayesian technique to derive estimates of 
the absolute magnitudes in U, B and V bands f o r each 
gala xy in our s a mple, following iKauffmann et aEI (|2003l ) 
and ISalim et al.! (|2005l ). Libraries of Monte Carlo realiza- 
tions of model star formation histories are generated be- 
tween < 2 < 0.5 in regular bins of Az — 0.001, using 
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the BC03 population synthesis code. Each hbrary contains 
25000 models with each star formation history being char- 
acterized with two components: an underlying continuous 
model with an exponentially declining star formation law 
and random bursts superimposed on this continuous model. 
The models also have metallicities and dust attenuation uni- 
formly distributed over wide ranges. The universal initial 
mass function (IMF) of iKroupal (|200ll ) is adopted. For each 
galaocy in our sample, we derive the U, B and V magni- 
tudes by comparing the observed ugriz SED to all the model 
SEDs in the closest redshift model library. The goodness 
of fit of each model determines the weight, oc exp(— 
which is assigned to that model when building the probabil- 
ity distribution functions (PDFs) of the restframe U, B and 
V magnitudes of the galaxy. We adopt the PDF-weighted 
mean values as our estimates of these quantities. Adopting 
the median of the PDF gives almost identical results for the 
autocorrelation function analysis. 



2.2 DEEP2 galaxy sample 

The DEEP2 Galaxy Redshift Surv ev jDavis et al.l 120031 ') 
utilizes the DEIMOS spectrograph ijPaber et al.l l2003l ) on 
the KECK II telescope. Targets for the s pectroscopic sam - 
ple were selected from BRI photometry (|Coil et al.ll2004al ) 
taken with the 12k x 8k mosaic camera on the Canada- 
France-Hawaii Telescope (CFHT). The images have a lim- 
iting magnitude of -Rab ~ 25.5. Since the _R-band provides 
the highest signal-to-noise ratio (S/N) among all the CFHT 
bands, the photometry in this band was used to select tar- 
gets for spectroscopic observation in the DEEP2. The CFHT 
imaging covers four widely-separated regions, with a total 
area of 3.5 deg'^. In fields 2 to 4, the spectroscopic sample 
is preselected using ( B — R) and (R — I) colors to eliminate 
objects with z < 0.7 (jPavis et al.ll2003 ^. Color and apparent 
magnitude cuts were also applied to objects in the first field, 
the Extended Groth Strip (ECS), but these were designed 
to downweight low redshift galaxi es rather than eliminate 
them entirely jWillmer et al.ll2006l '). 

We use data from the third data release of the DEEP2 
survejQ which contains spectra of about 50000 galaxies in 
the magnitude range 18.5 ^ -Rab ^ 24.1. The spectra have 
a resolution of ii ~ 5000. For this study we have selected a 
sample of 30546 galaxies from the DEEP2 DR3, with red- 
shifts of quality 3 or 4 and in the range 0.8 < z < 1.35. The 
median redshift of this sample is z = 1, which we use for 
comparisons with dark matter and model galaxy autocorre- 
lations. 

The derived galaxy parameters required in this work in- 
clude stellar mass (Af,) and restframe magnitudes in U and 
B bands. The procedure for estimating these parameters is 
exactly the same as the one above for the SDSS. In brief, 
the observed BRI SED of each galaxy in the DEEP2 is com- 
pared to a large grid of BC03 model SEDs, providing a max- 
imum likelihood estimate of the /-band mass-to-light ratio 
of the galaxy, as well as its restframe U and B magnitudes. 
A lKroupal (|200ll ) initial mass function (IMF) is adopted, as 
in the SDSS analysis. Our estimates of stellar mass and rest- 
frame magnitudes are statistically well consistent with those 
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from lBundv et all (|2006l ') and lWillmer et aL I (|2006D . Indeed 
we have repeated our clustering analysis using stellar mass 
and restframe magnitude estimates from these previous pa- 
pers, obtaining almost the same results. 

We don't consider the V band for DEEP2 galaxies, as 
at z ~ 1 this band is shifted well beyond the reddest band 
(the I band) that is observed. 

2.3 Semi-analytic model galaxy catalogues 

In this paper we compare our observationa l results to predic- 
tions from the galaxy formation model of lCuo et al.l (|201ll . 
hereafter Gil). This model was created by implementing 
semi-analytic models for baryonic astrophysics on merger 
trees encapsulating the e volution of the halo/subhalo popu- 
lation in the Millennium jSpringellbOOsI ) and Millennium-II 
iBovlan-Kolchin et al.ll2009l ') Simulation^ The Millennium 
was carried out in a cubic region 500 ft^^Mpc on a side 
with mass resolution ~ 10® M©, while the Millennium-II 
followed evolution in a region with 125 times smaller vol- 
ume, but at 125 times better mass resolution. The combi- 
nation of the two simulations allows galaxy formation to be 
studied over the full range of observed populations, from 
dwarf spheroidals to cD gala:xies. In comparison to earlier 
semi-analytic models from the Munich group, the treat- 
ments of supernova feedback, galaxy size, photoionisation 
suppression and environmental effects on satellite galaxies 
have been significantly updated, resulting in excellent fits 
not only to recent SDSS data on the luminosity and stellar 
mass functions of galaxies, but also to the recent determina- 
tions of the abundance of faint satellite galaxies around the 
Milky- Way. Most of the population properties for galaxies 
in the local Universe are reasonably well reproduced by the 
model, which also makes precise predictions for nonlinear 
clustering over the full observed range of scales, 10 kpc to 
10 Mpc. Comparisons between the clustering properties pre- 
dicted by the model and measurements of SDSS galaxy clus- 
tering show that agreement with SDSS data is quite good for 
masses above 6 x 10^" M0 and at separations above 2 Mpc, 
although the predicted clustering is up to 20% too high in 
some mass ranges. On smaller scales lower mass galaxies 
are predicted to be substantially more clustered than is ob- 
served. Gil suggest, but do not prove, that this may be a 
consequence of an overly large present-day fluctuation am- 
plitude (erg — 0.9) in the simulations. 



3 CLUSTERING MEASURES 

Following Papers I and II, we weight each galaxy in our 
samples to correct for incompleteness when computing our 
stellar mass (or luminosity) autocorrelations. The weights 
used for SDSS galaxies take into account three factors. The 
first is l//sp, where fsp is a spectroscopic completeness, de- 
fined as the fraction of the photometrically defined target 
galaxies for which usable spectra were obtained. The sec- 
ond weight, 1/Vij, is applied to each pair of galaxies, where 



^ Galaxy catalogues for this model and halo catalogues 
for the parent simulations are publicly available at 
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Vij = mm(ymax,i,Vmax,j) and Vmax,i IS the maximum vol- 
ume over which the ith galaxy would be included in the 
sample. This weight accounts for the fact that faint galaxies 
are not detected throughout the entire survey volume in a 
flux-limited survey. The final weight is the factor 1/ fcou.ij, 
which is also applied to galaxy pairs and appears in data- 
data counts only. This factor is a function of angular sepa- 
ration 6ij of the two galaxies and corrects for the effect of 
fibre collisions on small scales. 

Similarly, for each galaxy pair in the DEEP2 sa mple we 
also assign a weig ht following IWillmer et all (|2006l '): 



(1) 



where accounts for incompleteness resulting from the 
DEEP2 colour selection and redshift success rate. The sec- 
ond factor Vij is defined in the same way as above. Detailed 
de scription on the calc ulation of the weights can be found 
in IWillmer et all (|2006l V 

We have generated a random sample which has the same 
overall sky coverage and redshift distribution as the DEEP2 
sample. The spatial window function of the DEEP2 survey 
is applied to the random sample, which includes masking 
areas around bright stars and takes into account the vary- 
ing redshift completeness of the observ ed slit masks. The 
reader is re ferred to previous papers bv lCoil et al.l (|2004bl . 
I2OO6I . |2008| ) for detailed description of the window function. 
For each real galaxy we generate 100 sky positions at ran- 
dom within the entire survey region of the DEEP2 including 
all the four separate fields, and we assign to each of them 
the properties of the real galaxy, in particular, its values of 
K and Vmax- Since the four fields of the survey are widely 
separated, correlations in these properties in the r eal sample 
are wi ped out by randomizing in ang le. We follow I Coil etal] 
(|2006l ) to assign a redshift to each random point according 
to the redshift distribution averaged over all the fields in the 
data. 

The projected autocorrelation function of stellar mass 
or luminosity in a given band is then computed using the 
estimator described in Paper I (see their eqn. 2), in which 
the weights obtained above are applied. We estimate the 
autocorrelation functions using the entire real (or random) 
sample, for both SDSS and DEEP2. Errors on the mass 
autocorrelation in the SDSS are estimated from the scat- 
ter among the measurements from 20 mock galaxy cata- 
logues constructed from the Millennium Simulation using 
the same selection criteria as the real sample (see Paper I 
for details). The errors on the DEEP2 mass autocorrelation 
function come from the scatter among the measurements for 
the four separate fields of the survey. These errors should in- 
clude both the effect of counting noise and that of cosmic 
variance, which are impressively small in the SDSS due to 
the large size and volume of the sample. Following Paper 
I, we do not attempt to put independent error bars on the 
projected luminosity autocorrelations because, for a given 
survey, with our technique the set of galaxy pairs used to 
estimate each of these functions is exactly the same and as 
a result the noise fluctuations due to sampling and to large- 
scale structure are identical in all the luminosity and stellar 
mass autocorrelation estimates. 

It is important to point out that we have ignored the 
undersampling of DEEP2 galaxies on small scales due to 



the slit mask target-selection a lgorithm. Coil et al.l (|2006l ) 
use the mock galaxy catalogue o f lYan et al. I (|2004l ') to correct 
for the effect of this on their measurements of projected two- 
point autocorrelation functions, finding that the corrections 



are most significant on scales less than 



0.3/i-^Mpc. 



This scale is close to the lower limit {jp ~ 0.2/i~^Mpc) to 
which we plot our autocorrelation functions. We decided not 
to attempt to correct for slit mask effects, but one should 
keep in mind that our Wp(rp) measurements will be underes- 
timated fo r Tp < 0.3/i,~^Mp c, by roughly 20% on the smallest 
scales (see lCoil et"al]|2008l '). 



4 RESULTS 

4.1 Projected autocorrelation functions and bias 
factors 

In Figure [1] we show projected stellar mass autocorrelation 
functions Wp{rp), as well as projected luminosity autocor- 
relations Wp{rp) in the U, B, V bands. Results are plotted 
for the SDSS in the left-hand panel, and for the DEEP2 in 
the right-hand panel. In both panels Wp{rp) is shown using 
triangles with error bars, while the ■Wp{rp) are shown using 
lines. 

The left-hand panel of Figure[T]reproduces the results of 
Paper II where luminosity autocorrelation functions were es- 
timated for the five passbands of SDSS (ugriz). Luminosity 
is less strongly clustered than stellar mass on all scales and 
in all wavebands. Furthermore, the amplitude and slope of 
the luminosity autocorrelation function changes systemati- 
cally with wavelength, with the reddest band (the IZ-band) 
showing the highest amplitude and the steepest slope. As a 
result, y-band light clusters more similarly to stellar mass 
than the bluer bands plotted, but not as closely as the z- 
band light analysed in Paper II. This is expected because 
2-band light is known to be closely re lated to stellar mass 
indicators (e.g. iKauffmann et al.|[2003t ). Finally, all the au- 
tocorrelation estimates are well described by power laws. 
The best power-law fit to the stellar mass autocorrelations 
is shown in the figure, but in this paper we will not discuss 
such fits further, since they were a major topic of the two 
previous papers and the DEEP2 data do not allow autocor- 
relation shapes to be studied in as much detail as is possible 
at low redshift using SDSS. 

As can be seen from the right-hand panel of Figure [TJ 
the DEEP2 galaxies show very similar autocorrelation be- 
haviour to the SDSS. Although the measurement errors are 
relatively large due to the smaller sample size and the con- 
siderable depth of the DEEP2, all the systematic trends seen 
above for the SDSS hold also for the DEEP2. Furthermore, 
the slopes of the autocorrelation functions are consistent 
with remaining constant from 2 ~ 1 to the present, al- 
though their amplitudes have increased by a factor of about 
1.6 when, as here, they are measured in comoving coordi- 
nates. To be more quantitative, we have estimated the rela- 
tive bias factor as a function of projected separation between 
the SDSS and the DEEP2 by directly comparing their stel- 
lar mass autocorrelations. A linear fit indicates that the bias 
is well described by 6 = 1.6 — 0.09 log ]^Q(rp), with a reduced 
of 1.03. The rms scatter around the fit is 22.2% over the 
Tp range probed. Thus the slopes of stellar mass autocorre- 
lation from the two surveys are really very similar to each 
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Figure 1. Projected stellar mass autocorrelation functions in the SDSS/DR7 (left panel) and in the DEEP2/DR3 (right panel) are 
plotted as triangles with error bars and are compared to projected luminosity autocorrelation functions measured in the rest-frame U, 
B and V bands ( lines). Errors on the stellar mass autocorrelation function are estimated in the left panel from the scatter among 
measurements on 20 mock SDSS catalogues constructed from the Millennium Simulation, and in the right panel from the scatter between 
measurements on the four independent fields of the DEEP2 survey. The solid line in the left panel is a power-law fit to the SDSS data 
over the range 10h~^kpc < Tp < V)h~^ Mpc and is repeated in the right panel. It corresponds to a three-dimensional autocorrelation 
function ^*(r) = (r/ro)~^'** with tq = 6.1/i~^Mpc. The dotted line in the right panel is a power-law fit to the DEEP2 data over the 
full range probed, with index fixed to be -1.84, i.e. the same as what is determined for the SDSS data. The corresponding correlation 
length is tq = 4.3/i~^Mpc. 




fp [h ^Mpc] fp [h ^Mpc] 



Figure 2. Ratio of the projected stellar mass (triangles) and luminosity (lines) autocorrelation functions in the SDSS/DR7 (left panel) 
and in the DEEP2/DR3 (right panel) to the projected dark matter autocorrelation function at z = 0.07 (left panel) and at 2; = 1 (right 
panel) in the Millennium Simulation. Errors on the stellar mass autocorrelations are estimated as in Fig. [T] Errors on the simulated 
functions arc much smaller. 
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other. Our result is broadly consistent with previous studies 
of the evolution of galaxy two-point correlation functions us - 
ing data from the same surveys (e.g. ICoil et al.|[2006l . l20o"^ . 

These results are shown again in Figure [2] in terms 
of bias factors which we define as the ratio of the pro- 
jected stellar mass (or luminosity) autocorrelation to the 
projected dark matter density autocorrelation. The latter 
is obtained from the z = 0.065 and z = 1 snapshots of 
the Millennium Simulation and thus assumes the cosmolog- 
ical parameters of the simulation. As noted above, these 
redshifts are appropriate to characterise the mean depth 
of our SDSS and DEEP2 measurements. A maximum line- 
of-sight depth |A7r| — 40/i^^Mpc was adopted when com- 
puting the projected autocorrelation function for dark mat- 
ter in order to mimic our procedures for the real data. At 
Tp > 1.5/i~^Mpc all our estimates are consistent with bias 
being scale-independent. On these scales the overall bias at 
2 ~ 1 is higher by a factor of ~ 1.7 than at 2 ~ 0. This 
reflects the well-established result that large-scale structure 
has been growing more rapidly in the dark matter distribu- 
tion since z = 1 than in the galaxy distribution. 

At smaller separations the scale-dependence of the 
bias is strong, with a total range of a factor of 5 over 
10/i~^kpc< < l/i~^Mpc in the SDSS, and a factor of 
2 over 100/i"^kpc< < l/i^^Mpc in the DEEP2. Another 
obvious feature is a step at around 2 Mpc which is seen 
in all bias functions and at both redshifts. Considering the 
large error bars, we would say that the step feature in the 
DEEP2 bias function is significant only at about l.Scr (at 
0.5 h~^Mpc) to 2(T (at 0.7 /i~^Mpc) level. In order to quan- 
tify how closely the SDSS and DEEP2 bias factors approxi- 
mate each other, we have estimated the ratio of the two bias 
functions shown in Figure [2] as a function of projected sepa- 
ration. This can also be well described by a linear function 
6 = 0.58 + 0.08 logy)(rp), with a reduced of 0.4. The rms 
scatter around this line is 12.7%. 

Given the featureless Wp{rp) observed for stellar mass 
and light, the strong scale dependence of the bias factors 
on these scales is largely due to the much more pronounced 
features seen in the dark matter autocorrelations, mainly 
the remarkable change in slope at the transition between 
the one-halo term, where the pair counts are dominated by 
galaxy pairs in the same halo, and the two-halo term, where 
galaxy pairs are mostly in separate haloes, at a few Mpc (see 
figures 2 and 6 of Paper I for illustrations). As discussed in 
Paper I, this uncomfortable step feature suggests that the 
amplitude of mass fluctuations is too high in the Millennium 
Simulations, and that the shape of our bias functions might 
be used to estimate the value of the fluctuation amplitude 
parameter erg. We will come back to this point in the last 
section. 

4.2 Comparisons with a galaxy formation model 

In this section we compare the observed projected autocor- 
relations of stellar light and mass to predictions for these 
same statistics from the galaxy formation model of Gil, 
based on the Millennium and Millennium-II simulations. In 
their paper these authors already compared their model to 
SDSS clustering data, in particular, to the projected auto- 
correlations of stellar mass and of galaxies separated into 
passive and actively star-forming objects in a series of dis- 



joint ranges of total stellar mass. They found their model to 
reproduce observed clustering reasonably well, to a level bet- 
ter than about 20% at all separations for M, ^ 6 x 10^" 
and at rp > 2 Mpc for M. > 6 x 10^ M©, but to substantially 
overpredict the clustering of stellar mass at smaller separa- 
tions. Further comparisons between the SDSS and the model 
for autocorrelation functions in different intervals of stellar 
mass and optical colour showed this discrepancy to be due 
mainly to passive galaxies with stellar masses in the range 
6 X lO^'Af© < M* < 6 X 10^° Mq, indicating that the model 
predicts too many red, passive, satellites of this mass[f|. The 
authors suggested that this might indicate that the fluctu- 
ation amplitude erg = 0.9 adopted in the simulations is too 
high. 

Here we extend this study to 2 = 1 by comparing the 
model predictions with our DEEP2 results. The results are 
shown in Figure [3] (the right-hand panel), where we plot 
the ratios of observed to predicted Wp{rp) both for stellar 
mass (symbols with error bars) and for luminosity in the U , 
B and V bands (lines). For comparison, the corresponding 
results for SDSS are shown in the left-hand panel. At both 
redshifts the discrepancy between model and data is clearest 
in the stel l ar ma ss autocorrelations. As already found by 
IGuo et al] (|201ll ). the model Wp(rp) exceeds that observed 
in the SDSS on all scales, by only ~15% for rp > 2/i^^Mpc 
but by increasingly larger factors on smaller scales, reaching 
a factor of 3 at 10/i"^kpc. At redshift 2 1 our DEEP2 
results show a similar offset and are consistent with the same 
trend within their error bars. 

Remarkably, although the model predicts stellar mass 
autocorrelations that disagree with the SDSS data, its lu- 
minosity autocorrelations match observation much better, 
particularly in the rest-frame _B-band. Apparently the overly 
strong clustering of stellar mass is almost exactly compen- 
sated by the overabundance of passive satellite galaxies so 
that the light distribution is well reproduced. At z ~ 1 the 
predicted luminosity autocorrelations lie somewhat above 
the DEEP2 measurements in B but agree well at U . Given 
the error bars, the discrepancy at the longer wavelengths is 
only marginally significant. The fact that a physically based 
galaxy formation model can simultaneously agree with the 
clustering of light and disagree with that of stellar mass, yet 
agree with the abundance of galaxies as a function of both 
light and stellar mass (see Gil) demonstrates the complex- 
ity of the constraints on galaxy formation implied by precise 
observations of abundance and clustering. The agreements 
and disagreements with DEEP2 data at z 1 are very simi- 
lar to those with SDSS data at low redshift, suggesting that 
the model is treating the evolution and clustering of the 
galaxy population in a realistic way. 



5 SUMMARY AND DISCUSSION 

In this work we have estimated projected autocorrelation 
functions i«p(»'p) for the stellar mass of galaxies and for their 
light in the U, B and V bands, using data both from the 
third data release of the DEEP2 Galaxy Redshift Survey 

^ The "excess" of passive red satellite is well documented in the 
literature and seems to be a "common" problem for all recent 
galaxy formation models. 
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Figure 3. Ratio of the projected stellar mass (triangles) and luminosity (lines) autocorrelation functions in the SDSS/DR7 (left panel) 
a nd in the DEEP2/ DR3 (right panel) to predictions of these same quantities from the semi-analytic galaxy formation simulation (SAM) 
of lGuo et all 1I2OIII I at z = 0.07 (left panel) and at z = 1 (right panel). 
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Figure 4. Triangles with error bars indicate bias as a function of scale for stellar mass in the SDSS (left panel) and in the DEEP2 (right 
panel) assuming the Millennium Simulation cosmology but with varying values (as labelled) for the present-day fluctuation amplitude 
erg. For erg values other than 0.9, projected dark matter autocorrelations at 2: = 0.07 (for the SDSS) and z = I (for the DEEP2) are 
approximated by the functions found in the Millennium Simula tion at earlier times when the linear fluctuation amplitude matches that 
desired. Bias functions for stellar mass in the lGuo et al.l 1 I2OIII) galaxy formation model (which assumes erg = 0.9) at z = 0.07, 1 and 3 
are shown by solid lines, labelled by redshift. For clarity, error bars are only shown in the right-hand panel in the case erg = 0.8, and the 
triangles for each erg are connected by a dotted line. Only for erg = 0.8 or 0.7 is the bias function inferred from the data approximately 
as flat and featureless at the one-halo/two-halo transition as that measured in the model. 



and the final data release of the Sloan Digital Sky Survey 
(SDSS). We have compared these to projected autocorrela- 
tions of dark matter estimated from the Millennium Simu- 
lations at 2; = 1 and 0.07, the median redshifts of our galaxy 



samples, in order to investigate the bias of the clustering of 
stellar mass and light. 

The dependence of clustering on luminosity and colour 
for galaxi e s in the D EEP2 was previ ously studied by 
ICoil et al.1 (|2006l ) and ICoil et al.1 l|2008l) : the former in- 
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vestigates the luminosity-dependence of galaxy clustering, 
while the latter investigates the joint color and luminosity- 
dependence. These studies show that the dependence of clus- 
tering on color is stronger than on luminosity. These find- 
ings are qualitatively ve ry similar to those seen at low red- 
shift in S PSS data (e.g. IZehavi et al.ll2005l : iLi et al.ll2006al : 
IZehavi et al., ,201h ). Their conclusion was that the correla- 
tions between environment and galaxy properties which are 
well known at 2: = were, in fact, largely in place at z = 1. 
In this paper we have investigated these issues further by 
evaluating alternative two-point statistics, the projected au- 
tocorrelations of stellar mass and light, in a uniform way for 
the SDSS and the DEEP2, so that the two surveys can be 
compared directly to each other and to the predictions of a 
recent MS-based simulati on of the format ion and evolution 
of the galaxy population (|Guo et al.ll201ll ). These statistics 
have some advantages in that they allow the full observed 
samples to be used to characterize the distribution of stars 
at high precision over the comoving ranges 10 kpc to 30 Mpc 
in the SDSS, and 200 kpc to 10 Mpc in DEEP2. 

The most impressive result of this work is the remark- 
able similarity in behaviour between the SDSS and the 
DEEP2. The stellar mass and stellar light autocorrelations 
are, to within the still considerable error bars for DEEP2, 
identical in the two surveys except for an overall offset in 
amplitude by about a factor of two. Papers I and II have 
found that the autocorrelations of stellar mass (or lumi- 
nosity) are dominated by contributions from galaxies in a 
relatively narrow range in mass (or luminosity), around the 
characteristic value of the stellar mass (or luminosity) func- 
tion. In effect, this paper thus compares clustering of the 
population around M* (or L*) at the two redshifts. The 
similarity between the two surveys supports the conclusion 
of earlier analyses that the environmental dependence of 
galaxy properties seen in the local Un i verse was already in 
place at z = 1 (e.g. ICoil et al.l [2006L I2OO8I : ICooper et all 
I2OO6I : iMeneux et al] FoOd |2008| . [2009f ). This similarity, as 
well as the factor-of-two offset between the two redshifts, 
are well reproduced by the galaxy formation model of Gil 
which matches the clustering of stellar light very well and 
overpredicts the clustering of stellar mass to a similar extent 
at both epochs. 

The galaxy formation simulation fits the observed lu- 
minosity autocorrelations to within 30% on all scales, both 
at 2: = and z = 1, and it fits the stellar mass autocor- 
relations to about 15% at rp > 2 Mpc. On smaller scales 
the discrepancy in the stellar mass autocorrelations reaches 
a factor of 2 to 3. Correlations on large scales are pro- 
duced by galaxies in different haloes (predominantly the 
central galaxies of those haloes) and the tight correlation be- 
tween central galaxy properties and halo mass means that 
any model which fits the observed abundance of galaxies 
prod uces approximately correct correlations on these scales 
(e.g. lVale fc Ostrikerll2004l: IConrov et al.ll2006l: IWang et"al] 



l2006| - |More et al.ll2009l : lMoster et al.ll2010l : lGuo et al.l2010t ) 

Smaller scale correlations result from galaxy pairs residing in 
the same halo, so the discrepancy on these scales shows that 
the model has too many pairs of relatively massive galaxies 
which share a common halo. This may be explained by too 
large a value of cts in the Millennium Simulations, which re- 
sults in too many high-mass haloes which would then host 
these pairs. The much better agreement of the stellar light 



autocorrelations on these same scales, shows that the "ex- 
tra" satellite galaxies must have relatively little light for 
their mass, thus high stellar mass-to-light ratios. 

An interesting feature in the observed bias functions of 
Fig. [2] is the obvious step at the one-halo/two-halo transi- 
tion at about 2 Mpc which reflects the marked change in 
slope of the dark matter correlations at this point (see figs. 
2 and 6 of Paper I). We noted in Paper I that since physi- 
cally consistent models for the evolution of the dark matter 
and galaxy distributions show smooth bias behaviour over 
this separation range, this feature suggests that the ampli- 
tude of mass fluctuations is too high in the Millennium Sim- 
ulation, producing a one-halo/two-halo transition which is 
stronger in the simulated mass distribution than in the true 
mass distribution. If this interpretation is correct, the shape 
of our bias functions can be used to estimate the value of 
the fluctuation amplitude parameter erg. We illustrate this 
possibility in Figure |4] where we plot the bias determined 
from the clustering of stellar mass over the separation range 
200h-^kpc <rp< l^h'^ Mpc both for the SDSS {z = 0.07) 
and for the DEEP2 [z = 1). Here we use dark matter cor- 
relations measured in the Millennium Simulation not only 
at the median redshifts of the observed samples, but also at 
higher redshifts in order to represent (approximately) the 
correlations expected at the observed redshifts in cosmolo- 
gies similar to the Millennium cosmology but with lower 
erg. We choose earlier redshifts with linear fluctuation am- 
plitudes corresponding to ag, = 0.8, 0.7 and 0.6. The filled 
triangles show estimates of the stellar mass bias for the two 
surveys made in this way for four different values of erg. It is 
clear that as the fluctuation amplitude is decreased the step 
feature becomes less prominent in both surveys and the bias 
becomes a more strongly decreasing function of scale. 

Theoretical predictions for these stellar mass bias func- 
tions for z — 0.07, 1 and 3 are shown for the galaxy for- 
mation model of Gil as solid lines in Figure [l] As for the 
earlier galaxy formation models plotted in Paper I, none of 
these curves shows a feature at the one-halo/two-halo tran- 
sition, but their slope increases signiflcantly towards higher 
redshift. Comparing with the observational results, it is clear 
that the the best agreement in amplitude occurs for ag = 0.8 
whereas the overall shape agrees best for somewhat lower 
amplitudes. It is remarkable that these results hold both 
for the SDSS at z = 0.07 and for the DEEP2 at z = 1. 
Hence the models appear to describe the overall evolution 
of the stellar mass autocorrelations surprisingly well, once 
the cosmology is corrected to a lower erg value, similar to 
that indicated b y most other recent analyse s of CMB and 
other data ( e.g. Komatsu et al.l |2009| . [2OI0I ) (However see 
IWang et al. I l|200a )" who argued that lowering ag would not 
always solve the problem.). 

In this paper we have only considered the diagonal er- 
rors on our autocorrelation function measurements, since we 
wish only to characterise the current situation, presenting 
a side-by-side comparison between the SDSS and DEEP2 
in terms of the observed autocorrelations, the bias relative 
to dark matter, and a comparison with the current semi- 
analytic model. To go beyond this and to make rigorous 
model fits to the full ensemble of Wp{rp) estimates would be 
a major undertaking that is not possible here. 
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